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Abstract We establish a new self-consistent system of equations accounting for a non- 
minimal interaction of gravitational, electromagnetic and axion fields. The procedure is based on 
a non-minimal extension of the standard Einstein-Maxwell-axion action. The general properties of 
a ten-parameter family of non-minimal linear models are discussed. We apply this theory to the 
models with pp-wave symmetry and consider propagation of electromagnetic waves non-minimally 
coupled to the gravitational and axion fields. We focus on exact solutions of electro dynamic equa- 
tions, which describe quasi-minimal and non-minimal optical activity induced by the axion field. 
We also discuss empirical constraints on coupling parameters from astrophysical birefringence and 
polarization rotation observations. 
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1 Introduction 



The topic of pseudoscalar-photon interaction and axion theory is a very attractive branch 
of modern physics, and there are two reasons for this interest. First of all, it is an interest 
to axions as hypothetical particles, which appear in the context of the strong CP-violation 
problem and spontaneous breaking of symmetry in the early universe (see, e.g., [TJ [2] for 
review and basic references) , and which are considered as one of the Dark Matter candidates 
(see, e.g., (3j H] and references therein). The direct detection of the axions is formulated to 
be one of goals of modern experiments in High Energy Physics (see, e.g., [5]). 

The second reason is connected with an interest to the theoretical concept of pseudoscalar 
field cf) (axion field) associated with pseudoscalar-photon interaction. This occurs first in the 
study of electromagnetism and equivalence principles [6J [8] . This concept is a base for the 
so-called Einstein-Maxwell-axion theory, Einstein-Maxwell-dilaton-axion theory and their 
generalizations (see, e.g., [6]-[T8]). 

One of the most significant item in the theory of axion fields is the problem of interac- 
tion of electromagnetic and pseudoscalar (axion) fields, the well-known application of the 
photon-axion coupling being the effect of polarization rotation [6j [191 120] in the observa- 
tions of Cosmic Microwave Background (CMB) radiation (see, e.g., [21] [22] and references 
therein). Several groups are working on the experiments with vacuum birefringence and 
vacuum dichroism to look for the photon-axion coupling [23]-[26j. 

It is well-known that the standard effect of optical activity induced by the axion field 
takes place when the pseudoscalar field has a non- vanishing four-gradient [8j [27], i.e. when 
dk4> 7^ 0. This effect can be described in the framework of classical electrodynamics in terms 
of linear constitutive equations [Hi E7J [28] , linking an excitation tensor H lk and the Maxwell 
tensor F mn . When the pseudo-scalar field is constant, <fi = 0o, and thus does not contribute 
to the vacuum Maxwell equations, new possibilities for the axion-photon interaction exist: 
one of them is discussed in the work [29], which focuses on reflection and transmission of 
a wave at an interface between two media. We discuss here another possibility related to 
non-minimal optical activity induced by photon-axion coupling, which also can occur in the 
case of a constant pseudoscalar field. 

The standard description of the non-minimal coupling of the gravitational field with 
scalar, electromagnetic and gauge fields is based on the introduction of specific cross-terms 
into the Lagrangian, which contain the Riemann tensor Rtkmn, the Ricci tensor, Rkn, and 
Ricci scalar, R, on the one hand, and the corresponding fields and their derivatives, on the 
other hand. The theory of non-minimal coupling is elaborated in detail for the scalar field 
(real and complex fields, as well as Higgs multiplets) (see, e.g., [30] for review and references). 
Special attention in these investigations is focused on two models, the first of them has the 
£i?$ 2 coupling, and the second has the so-called non- minimal derivative coupling |31j-[33j. 
The study of a non-minimal coupling of gravity with electromagnetic field started in [34] and 
this theory has been developed by many authors (see, e.g., [3S]-|1H])- The generalization of 
the idea of non-minimal interactions to the case of torsion coupled to the electromagnetic 
field has been made in [49, 50] . Exact solutions in the framework of non-minimal Einstein- 
Yang-Mills and Einstein- Yang-Mills-Higgs theories are discussed in |51j-[56j. 

In this paper we focus on the non-minimal generalization of the Einstein-Maxwell-axion 
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theory. The derivation of master equations is fulfilled in analogy with non-minimal scalar 
field theory and non-minimal electromagnetic field theory. The items which are specific for 
the pseudoscalar field theory are discussed in detail. The paper is organized as follows. In 
section 2, we discuss a formalism of non-minimal Einstein-Maxwell-axion theory; a brief 
introduction to the minimal theory is given in subsection 2.1; the non-minimal extension of 
the Lagrangian and the derivation of the corresponding master equations are presented in 
subsection 2.2; subsection 2.3 contains the description of the model in terms of constitutive 
equations. In section 3 we apply the obtained master equations to the models with pp-wave 
symmetry and discuss examples of the exact solutions for the gravitational, pseudoscalar 
and electromagnetic fields. In section 4 we consider the propagation of test electromagnetic 
waves coupled to pseudoscalar field in the pp-wave background and focus on the effect of 
non-minimal optical activity induced by the axion field. In section 5 we summarize the 
results. 



2 Non-minimal coupling of gravitational, 
electromagnetic and axion fields 

2.1 Minimal model as a starting point 

In order to explain the novelty of our approach, let us first introduce the case of gravitational- 
electromagnetic-axion fields minimally coupled. The simplest minimal action functional is 



S(M)= / d xy/^g 



^ + I F - Fmn+ l0 F *- Fmn _^-v m 0V ri 0+m2 ) 2 
K 2 2 y ' 



A) 



where g is the determinant of the metric tensor g^, V m is a covariant derivative, R is the 
Ricci scalar. The Maxwell tensor F mn is given by 



Fmn — V ' m,A n V ' n A m , (2) 

where A m is an electromagnetic potential four-vector; F* mn = \e mnvq F m is the tensor dual 
to F pq ] e mnpq = -j=F mnpq is the Levi-Civita tensor, E mnpq is the absolutely antisymmetric 
Levi-Civita symbol with E 0123 = 1. The dual Maxwell tensor satisfies the condition 

\/ k F* ik = . (3) 



The first term in the brackets is the Hilbert-Einstein Lagrangian; the second term is the 
standard Lagrangian for an electromagnetic field; the third term is the pseudoscalar-photon 
interaction Lagrangian [6j [7J [8] ; the fourth and fifth terms constitute the pseudoscalar La- 
grangian. 

The symbol <fi stands for a pseudoscalar field; this quantity is dimensionless providing the 
terms \F mn F mn and \(j)F* mn F mn to have the same dimensionality. The axion field itself, $, 
is considered to be proportional to this quantity $ = ^o<p with a constant \l/o- In order to 
simplify the calculations we replace the potential of electromagnetic field A\ with the term 
Ai/^o, and consider the constant k to be equal to k — S7T( ^o _ term mtx) is proportional 



3 



to a (hypothetical) mass of an axion, m(A) = 27rc m( ax i on )//i; h is the Planck constant. We 

use the signature H . 

The variation of the action functional (Op) with respect to the four-vector potential Ai 
gives the minimal (vacuum) Maxwell equations 



F ik + <f)F 



*ik 







which can be transformed into 

V k F ik + F* ik V k (j) = 
due to the equation (j3J) . Equivalently, equation (jlj) can be written as 

V k H ik = , # i/c = C ifcmn F m „ , 



where C mn is the constitutive tensor 



C 



ikmn 



im kn in km\ , / ikmn 

9 9 -9 9 J + 



(4) 



(5) 



(6) 



(7) 



and H tk is the excitation tensor [HI [27] . Minimal equations for the axion field can be obtained 
from the action ([T]) by the variation with respect to the pseudoscalar field </>, yielding 



W k W k +m 2 (K) 



1 

T 



F*mn rp 



Minimal equations for the gravitational field obtained by the variation of ([TJ) with respect 



to metric g lk are 



Rik — -R9ik — K 



T (EM) T (A) 
1 ik ^ 1 ik 



Here T^ U ' is the standard stress-energy tensor of pure electromagnetic field: 



T: 



(EM) 



■ik 



-auF mn F — F F, m 

^ilik ± ± mn ± ira 1 fc j 



(9) 



(10) 



and given by 



(11) 



is the stress-energy tensor of the pseudoscalar (axion) field 0. 



2.2 Non-minimal extension of the model 
2.2.1 Lagrangian and non- minimal susceptibilities 

We consider the total action functional as a sum of minimal and non-minimal contributions: 

S = S(m) + S(NM) ) (12) 
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where S(m) is given by (0Q). The non-minimal contribution <9(nm) is, generally, a nonlinear 
function of all independent invariants containing Riemann tensor Rikim and the dual ones 
*Rikim, Rikirrn Ricci tensor R mn and Ricci scalar R in convolutions with the Maxwell tensor 
and its dual, as well as with derivatives of pseudoscalar field. The set of such invariants for 
pure electromagnetic field was discussed in [17]; then it was supplemented by the invariants 
containing scalar Higgs fields in [531 E31 EE] . For the axion field, this set of invariants can be 
constructed analogously, and here we focus only on the following Lagrangian linear in the 
curvature: 

Scnm)= J d A x^g~ \^ kmn F lk F mn + l -x l t^ n (t> F, fc F,; n -^V m 0V n 0+r/ { A)i?0 2 } • (13) 

The quantity TZ lkmn is a non-minimal three-parameter susceptibility tensor [17], which has 
a form 

n ikmn = qiRg ikmn + ^kmn + q ^ R ikmn ? ^ 

where 

g ikmn _ _^ g im g kn_ g in g km^ ? ( 15 ) 

U ikmn = l( R im g kn_ R in g km +R kn g im_ R km g inj _ ( lg ) 

The constants gi, g2 and are non-minimal parameters describing the linear coupling of 
the Maxwell tensor F mn with curvature [37] ■ The quantity X(A) n given by 

x f^ n =Q 1 Rg lkmn +Q 2 W kmn +Q3R ikmn , (17) 

where Qi, Q2 and Q3 are also constants, is the non-minimal susceptibility tensor describing 
the linear coupling of the dual tensor F^ n with curvature. As in the previous case, the 
combination 4>F^ n gives the tensor quantity. The tensor 

W$ = + F nl R m ^j + r] 2 Rg mn + r] 3 R mn (18) 

describes a non-minimal susceptibility for the pseudoscalar field in analogy with [53], but in 
this case it contains a term linear in the Maxwell tensor. This term describes effects analogous 
to the so-called derivative coupling in the non-minimal scalar field theory [311 [32j [33] . As for 
the tenth coupling constant 77(A) , it is a direct analog of the well-known coupling constant £ 
in the non- minimal scalar field theory (see, e.g., [30] for review and references). 

The tensors TZ lkmn and X(aT"> defined by (1141) and (1171) . are skew-symmetric with respect 
to transposition of the indices i and k, as well as m and n. In addition, the following relations 
take place 

sTjikmn srymnik ikmn mnik /i n\ 

><■ ~ K > X(A) - X(A) > l iy J 

which guarantee that the model under consideration does not contain the solutions of the 
skewon type [27J. The tensor 9ft^™ is explicitly symmetric, i.e. = 9Rya). 
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2.2.2 Non-minimal electrodynamic equations 

Electrodynamic equations, which correspond to the Lagrangian (Tl2"i) with ([[]) and ( !T3l) . are 
linear and have the standard form 

W k H ik = P . (20) 

The excitation tensor H lk and the Maxwell tensor F mn are linked by the linear constitutive 
law 

H ik = F lk + n ikmn F mn + [0 (F* ik + xfK) nF *mn)] • (21) 

The second term on the right-hand side is the curvature-induced polarization-magnetization, 
appearing in the non-minimally extended pure Einstein-Maxwell model [47J. The contribu- 
tion detailed in square brackets is the one of the axion type, the terms and F* lk enter the 
equations in the multiplicative form only (see, e.g., J8j|35] for more details). The non-minimal 
axion contribution 4>X{A) n F^ in is a new term proposed here. 
The four-vector of an effective electric current 

P = i^V fc [(i? fcm V*0 - R im V k <p) V m 0] (22) 

is proportional to the coupling parameter rji and is due to the first term in ([TBI) . This 
four-vector does not contain Maxwell tensor and satisfies the conservation law VjP = 0. 
It contains both the Ricci tensor and the pseudoscalar field; thus, it describes the electric 
current induced by curvature and axion field derivatives. 

2.2.3 Non-minimal equation for the pseudoscalar field 

Non-minimally extended master equation for the pseudoscalar takes the form 

V m [(g mn + RJS) V n 0] + [mf A) + V{A) R\ = --F^F^ - \ X f^ n F tk F* mn , (23) 

where 9ft^ and x\ k A) n are gi ven by f|T8|) and f|T7|) . respectively. This equation is a non-minimal 
generalization of (JS]). 

2.2.4 Non-minimal equations for the gravitational field 

Variation of the action functional f lT2|) with ([1]) and f|T3|) with respect to g tk gives the non- 
minimally extended equations for the gravitational field 

(fl»~«9.) d+«e) =* [rr ) +^ A '+rr EM, +^ NMA> ] • (24) 

Here the scalar stands for the quantity 

ee 77(A)0 2 +| qiF mn F mn + l - Q 1( j ) F^ n F mn + (^3-772) V m 0V"> , (25) 

with the tensors T^ M ^ and given by ( ITUl) and respectively. The non- minimal 

extension of the stress-energy tensor contains two contributions: first, 7 1 .( fc NMEM ) describing 
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pure non-minimal electromagnetic part (see, e.g., |U] for details); second, the non-minimal 

'ik 



axion part 7"[ NM . These tensors can be specified as follows: 



r' 



(NMA) 



T (NMEM) 
1 ik 



-(2) 



-(3) 



(4) 



-(5) 



-(6) 



ik 



where 



-(7) 
ik ' 



T, 



(1) 



ik 



RFi m F k m + 



1 r 



,(2) 



- ik 



ViVfc-^V'V, [F mn F r > 



V m V, (F mn F l n ) - R lm F mn F l % 



(26) 
(27) 

(28) 



-F n (RuF kn + R kl F m ) - R mn F m F kn - - V m V m {F m F k < 
+ l -V l [v< (F kn F ln ) + V fe (F m F ln 

Tift) 1 rymnls rp rp rpls / p nn 1 TP n p \ 

1 ik — ~^9ikti r mn ri s — -r n kn i s + r k tUnls) 



1 



V m V n [F. + F fc »F 4 



'ifc 



ViV fc -^V'V,J [0i^ n F m "] --g lk R<pF* mn F r 



'ik 



i4*F (RimF kn ~\~Rkm,Fi n ) ~\~ 



+ { V, [0 (F fcn F^+F fc * n F ta )] +V fc [0 (F in F* ln +F* n F ln )\ } 

- J^V n V m [0 (Fnl F *r ni+F *nl F m i 



-\v m V m [<f>(F in F™+F* n F k ")} 



4 4) = -g* (Rn-V l V n ) (F" m V m 0V,0) +^V/0 (F, n V k </>+F k n Vi<f>) + 



(29) 

(30) 
(31) 



(32) 
(33) 



+\ V'V, [V m (F™ V fe 0+F™ V^)] +^V ? [V, (F fc m V m 0V,0) +V fc (F i m V m 0V^)] + 



+ 4 V m [V, (F mn V k <pV n <P) +V k (F m "V,0V n 0)] + 

+^F mn (R m V k (f ) +R kn V l <f ) ) V m 0+^ {RTFl+R^F^ V m 0V n c 
TTf = FV,0V fc + (^V„V™ - V,V fe ) [V m 0V m 0] , 

4 6) = v m [i?rv fc + ^rv^] - i^ fc v m 0v">+ 



(34) 
(35) 
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+ i<7i*V ro V„ [V m 0V"0] - V m [V m (V, V fc 0)] . 

4 7) = (V,V fc -^V m V m )0 2 . 
Straightforward calculations show that the following identity takes place: 



v fc {(i+ K e)- 1 



T (EM) T (A) T (NMEM) r (NMA)l 1 _ n 
1 ik ' J ik ' 2 ik <~'ik J J — 



(36) 
(37) 

(38) 



i.e. the total effective stress-energy tensor written in the braces is a conserved quantity. 



2.3 Non-minimal constitutive equations for electromagnetic field 

Relation ( )2T|) is a linear constitutive equation [57J |27J [8] of the following type: 



H 



ik 



C 



ikmn : 



where 



Qikmn giKmn _|_ ^^gWcWMl _|_ 

Here we use the standard definitions for the right and left dualization 



j ikmn 



*ikmn , * ikmn 
A. (A) ' A. (A) 



*ikmn 
X(A) 



i*(A) ' 



-pq 



* ikmn 
X(A) 



ifc pqmn 
■ e pqX(A) 



(39) 



(40) 



(41) 



The tensor C lkmn describes the linear response of the material to the electromagnetic field 
action and contains the information about non-minimal dielectric permittivity and magnetic 
impermeability, as well as about the non-minimal magneto-electric coefficients in analogy 
with the standard continuum electrodynamics [57] [58 | [59]. 

Let us mention that in many works (e.g., in [27], |60j-[62j), a constitutive tensor density 
j^ikmn _ ^j—g Qikmn j g usec } f or the description of linear response instead of true tensor 
C lkmn . In terms of the quantity j^ kmn the condition for no birefringence (no splitting, no 
retardation) for electromagnetic wave propagation in all directions in the weak field limit 
gives ten constraint equations on the components of constitutive tensor density (see, e.g., 
[60j-[63]). With these ten constraints ■^ tkmn can be rewritten in the following form 



X 



•-ikmn 



km 



i „ rpikmn 



(42) 



where H = det(Hik) is a determinant of an effective metric H^, which generates the light 
cone for electromagnetic wave propagation; ip is some dilation factor; (p differs from our <p 
by the coefficient two <p — 2<P- I n case when (j4"2l) is not satisfied, birefringence will occur. In 
case when ( )42|) is satisfied, a pseudoscalar field ip will give a polarization rotation (an optical 
activity induced by ip). The effects of birefringence and polarization rotation described in 
terms of constitutive tensor density are studied in the 1970's and 1980's. In particular, 
constraints on birefringence from pulsar signal observations give the following estimates for 
the birefringence part Sx lkmn of j^ kmn in a weak field: 



5x 



ikmn 



< 10 



-14 



10 



-16 



(43) 
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Using the medium velocity four- vector U" 1 , normalized by U l Ui = 1, one can decompose 
Qikmn un iq Ue iy as 

Qikmn 1 ^imjjkjjn ^injjkjjm _|_ ^.knjjijjm ^.kmjjijjri^ _j_ 

■V iH {^~ l )isV mns + V iM ( um Vi n ~ V n v x m ) + r7' mn (C/V 2 fc - Ef*i/,*)1 . (44) 



Here e im is the dielectric permittivity tensor, (// l ) pq is the magnetic impermeability tensor, 
and v p m is the tensor of magneto-electric coefficients. These quantities are defined as follows: 

£ 2C UfJJ n , (/! )pg ~^1]pikC Vrnnq i 



" - r ha fif kmn U n = f4C mfc % np . (45) 
We use the symbols r\ mn i and A lh for the tensors 

Vmm = e mnls U s , A tk = g lk -U l U k , (46) 

orthogonal to C/\ The tensors and (/i" 1 )^ are symmetric, but z/^ is in general non- 
symmetric. These three tensors are orthogonal to U\ 

e lk U k = 0, {^ 1 ) ik U k = 0, v l k U l = = v, k U k . (47) 

Using expression fj40|) one can calculate the tensors e im , (/i _1 )j m and z/ pm explicitly. The 
symmetric dielectric permittivity tensor 

£ im = A im [l +qiR+q2R P* UpUq ] + g 2j R w A lp A m H 
+2q 3 R ipmq U p U q +Q 3 (f) U p U q (*RiPmq + R*i^ ? ( 48 ) 

as well as the symmetric magnetic impermeability tensor 

(jjT 1 )*" = A™ [l+ gi R+q 2 R pq U p U q ] -q 2 R pq A^A mq - 

-2q 3 *R* ipmq U p U q +Q 3 <f) U p U q (* R ^i +R »p^ (49) 
do not contain coupling parameters Q\ and Q 2 . The cross-tensor 



q 2 V pmk RkiU l + 2q 3 *R plmn U l U n + 



A 



pm 



Qi+~Q 2 } B 



+Q 3 <pUiU n (*R*Plmn_ RP lmn\ ^ ^ 



which gives the optical activity effects (see, e.g., [59]), contains two explicitly distinguished 
parts. The first part does not contain pseudoscalar field 0, is linear in the parameters q 2 and 
#3 and describes optical activity induced by direct non-minimal interaction between electro- 
magnetic and gravitational fields (such a model was analyzed in [65, 54J). It is important to 
mention that this part of the cross-tensor v pm contains both: skew-symmetric and symmetric 
terms (see the terms with r] pmk and * R plmn UiU n ) . The second contribution to u pm , which is 
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proportional to the pseudoscalar field 0, is symmetric with respect to transposition of indices 
(see the terms, containing the projector A pm and the tensor *R*P lmn —j{P lmn ^ it describes 
an optical activity induced by the axion field, the term A pm (j) relates to the minimal effects, 
the terms proportional to Qi, Q2 and Q3 relates to non-minimal effects. It is well-known 
that the non- vanishing cross-tensor u lk indicates that the medium is optically active, and the 
rotation of the Faraday type takes place in the course of electromagnetic wave propagation. 
Thus, one can see directly from fl50|) . that the non- minimal interaction between electromag- 
netic and axion fields produces curvature-induced optical activity of a new type. Below we 
consider this effect explicitly by the example of model with pp-wave symmetry. 
In table 1, we summarize briefly the information about non-minimal coupling constants. 





Term in the 
Lagrangian 


Physical meaning 


Qi 


1 td pmn p 

2 ± mn 


NMEM susceptibility linear in the Ricci scalar 


I2 


TDmn p p k 


NMEM susceptibility linear in the Ricci tensor 




1 TDikmn p p 

2 n, 1 ik-T mn 


NMEM susceptibility linear in the Riemann tensor 


Qi 


^RF™F* mn 


NMEM susceptibility induced by the axion with the Ricci scalar 


Q2 


<\>R mn F mk F* n k 


NMEM susceptibility induced by the axion with the Ricci tensor 


Qs 


^(f)R lkmn F ik F^ n 


NMEM susceptibility induced by the axion with the Riemann tensor 


Vi 


-R^F^VmCjNnCj) 


NMEM current induced by the axion field gradient 




-RV m (j)V m (f) 


NM axion-graviton derivative coupling with the Ricci scalar 


V3 


-R mn V m <j)V n <j) 


NM axion-graviton derivative coupling with the Ricci tensor 




Rcf> 2 


NM correction to the mass square of the axion 



Table I. Ten non-minimal (NM) coupling parameters are divided into four subgroups: the first 
Qi-, Q2, Q3', second Qi, Q2, Q3] third ^1,772,%; and fourth ?7(a)- In the second column the terms in the 
Lagrangian are given in front of which the corresponding coupling parameters are introduced; the 
parameters of the first subgroup introduce the terms without the pseudoscalar field (ft; the param- 
eters of the second subgroup relate to the terms linear in <j); the terms indicated by ^1,772,773, are 
quadratic in the four-gradient of (p; and finally, t/(a) introduces the term quadratic in (p. In the last 
column, we point out the physical meaning of these non-minimal terms, based on decompositions 
of the constitutive tensors for the electromagnetic (EM) and pseudoscalar fields. 



3 Non-minimal models with pp-wave symmetry 
3.1 Reduced non-minimal field equations 

We consider the line element for the spacetime with pp-wave symmetry in the form [66] 

ds 2 = 2dudv - L 2 {cosh 2 7 [e 2/3 (dx 2 ) 2 + e" 2/3 (dx 3 ) 2 ] + 2 sinh 2 7 dx 2 dx 3 } , (51) 

where w= ct ^| 1 and v= ct ^ are the retarded and the advanced times, respectively; /3(u) and 
7(ii) are functions of retarded time only, they determine the polarization of gravitational 
wave field. The case when the gravitational pp-wave field is characterized by one polarization 
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only (7 = 0) is investigated in detail in [M]. Various properties of pp- waves are discussed, 
e.g., in J67J E3 EH]. For the metric (IBTj) . the Riemann tensor has only three non- vanishing 
components R 2u 2u, R3u3u and R 2u 3u- The Ricci tensor has only one (generally) non- vanishing 
component R uu and the Ricci scalar is equal to zero identically R = 0. For this metric the 
three Killing vectors, which form an Abelian subgroup of the group G5 [69J, are 

C(v) = C(2) = $2, C(3) = 5 3 • ( 52 ) 

The Killing vector £Ls is a covariant constant null vector orthogonal to an d C(3)- We sug- 
gest here that the axion and electromagnetic field potentials inherit the spacetime symmetry, 
i.e. the Lie derivatives satisfy the following relations: 



£ 











(53) 



where (a) take the values (v), (2) and (3). As a consequence of ( )53l) . the pseudoscalar field and 
electromagnetic field potential depend on the retarded time only, = <f>(u), A m = A m {u). 
Taking into account the Lorentz gauge VkA k = 0, we can choose the potential in the form 



A u = A v = 0, A 2 = A 2 {u) , A 3 = A 3 (u) 



(54) 



The next step is to verify that pp-wave symmetry is admitted by the total system of non- 
minimal master equations. Direct calculations show that for such an ansatz about symmetry 
the non- minimal equation for pseudoscalar field 4>(u), i.e. the equations ( 1231) with (fl8l) . are 
satisfied identically for arbitrary Qi,...r] 3 and r/(A), when the axion field is massless, i.e. 
m (A) — 0. Also, the non-minimal electrodynamic equations f )20|) with (|2ip and ( 122]) are 
satisfied identically for arbitrary coupling constants, when the potential takes the form ()54p . 
Finally, the non- minimal equations for the gravitational field (j24j) with f l25"l) - fl3"Tj) can be 
reduced to one equation only, namely 



L" 



+ (f3'f cosh 2 2 7 + (7 



/\2 



1+^(A)0 2 ] = (0') 2 (l+2?7(A)) +2r/ (A 



+ 



+ 



L 2 



|cosh2 7 (A! 2 e~P) 2 + (A^y -2sinh2 7 A' 2 A' 3 j 



(55) 



The prime denotes here and below the derivative with respect to the retarded time u. The 
right-hand side of this equation includes, as usual, the minimal contributions from the tensors 



T ifc and T ifc . Only two non-minimal contributions came from the term l~ k f l37j) and from 
the term k7/(a)</> 2 in the expression for O fl25l) . The non- minimal coupling constants q±, q 2 , 
<?3, Qi Q2, Q3, V2, V3 are non- vanishing, but they happen to be hidden in this pp-wave 
model. Thus, in the presented model, six functions (p(u), A 2 (u), Az(u), (3(u), j(u) and 
L(u) are linked by one equation ( 151)1) only, and this model admits very wide possibilities in 
searching for exact solutions. 
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3.2 Exact solutions of the pp-wave type: 
quasi-minimal models with 77(A) = 

When the non-minimal parameter 77(A) is vanishing, we deal with a model which looks like 
minimal, i.e. the key equation for the pp-wave gravitational field 



2L 2 



L" 
cosh 27 



(/3') 2 cosh^2 7 + { 1 'Y+-( ( j ) 'Y + 



A2 



K 



L'\2 



G4< 



-2sinh2 7 A' 2 A' 3 



(56) 



does not contain non-minimal parameters at all. That is why we indicate this case with 
arbitrary but hidden parameters qi, 92, (73, Qi Q2, Q3, ^1, ^2, V3 an d vanishing r/(A) as quasi- 
minimal. Clearly, the quantity ( — x) * n ^ ne left-hand side of (l5"6"j) has to be non- negative. 
This property of L(u), which has initial value L(0) = 1, provides L(u) to vanish at some 
moment u = u* {L(u*) = 0)(see, e.g., [68]). This moment of the retarded time relates to 
the well-known singularity, the physical sense of which is discussed in detail in many books. 
In [BS], a number of papers are quoted, in which exact solutions of (|56j) are obtained for 
the pure electromagnetic source (0 = 0). In addition, exact solutions of the model of this 
type are known for vanishing Ak, when the (pure)scalar field $ is the source of the pp-wave 
gravity field (see, e.g., [TOJ ) . Finally, when Ak = and <fi = 0, we deal with the so-called 
pure gravitational pp-wave (see, e.g., [Ml EU |68j [69] ) . We would like to list here only three 
examples from the wide collection of minimal solutions of the equation ( 15 6 p . which relates 
to the models with symmetric spacetime of the pp-wave type. The models with symmetric 
spacetime are characterized by the condition ([67, 69J) 







In particular, when j(u) = 0, one obtains from ( 15 7p that 



R 2 — 
rl u2u — 


L" 

T 


+ 






= \i 


p3 _ 
rt u3u ~ 


L" 

T 




+ 


2^-/3' + f3" 


= A 2 



where Ai and A2 are some constants. The summation of ( l58l) and ( l59l) yields 

Ai + A 2 . 



T + (/3 ' )2 



(57) 



(5? 



(59) 



(60) 



Thus, for the symmetric spacetime, the key equation for the (quasi) minimal model can be 
reduced to 



-(Ai+A 2 ) = 

K 

In more appropriate terms, when 



(A'e^) +{A^y 



(61) 



F(u)G(u), 2/3 = log - 



(62) 
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(158]) . (!59|) give two independent equations 

F"(m) + A 1 F = 0, G"(m) + A2G = 0, (63) 
which can be easily solved if the signs of the constants Ai and A 2 are fixed. 

(i) First example 

The first model of such type relates to the pure gravitational wave, i.e. to the case when 
Ruu — 0- Suggesting that Ai = — A 2 = v 2 , we obtain immediately from ( 163]) the well-known 
Petrov's solution 



ds 2 = 2&U&V— cos 2 vu(dx 2 ) 2 — cosh 2 vu(dx 3 ) 2 , (64) 

with 

F{u) = cos vu , G{u) = cosh vu . (65) 

This solution is admissible, when 0=0o=const, Ak=0. The quantity y/—g = cos vu cosh vu 
vanishes at u* = n/2v. 

(ii) Second example 

When A2 = A3 = 0, (3 = 7 = and function <p(u) is linear in the retarded time, i.e. 
(f)(u) = 0o + wm, then ( 156]) gives 

L{u) = cos y^wM , L(0) = 1, L'(0) = 0. (66) 
For such a metric one obtains that 

R 2 u2u = r3 u3u = — j- = -y- = Ai = A 2 . (67) 
The first zero of the function L is w = u* 



w\/2k 

(iii) Third example 

When Ai = v 2 , Ai = — fi 2 and = 0, there is an exact solution 

F(u) = cosvu, G{u) = cosh/iu, (68) 

with 

v 2 = v 2 + ^. (69) 
The electromagnetic field in this model is presented by the following potentials: 

F/Q 

A2(u)=—— — [u sin uu cosvu— v sin vu cos uu] , (70) 

\uj l —v l ) 



Az{u)= ° [u sin oju sinh fxu — to cos cou cosh fiu\ . (71) 



Eg 

This electromagnetic wave can be considered as a circularly polarized wave, since the physical 
components of the electric field satisfy the following conditions: 

E\ = —(E 2 E 2 + E 3 E 3 ) = -g^A' a A' p = E 2 , (a, /? = 2, 3) . (72) 

When 77(A) 7^ 0, the quantity ( — "x) can change the sign, since the term r/(A)(0 2 )" in the right- 
hand side of ( 155]) is not positively defined in general case. This means that the singularity 
L = can, in principle, be avoided. 
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3.3 Exact solutions of the pp-wave type: 
Non-minimal models with 77(A) 7^ 

Consider now exactly solvable models with non- vanishing coupling constant 77(A)- 
3.3.1 Regular model 

Let us consider a model with L(u) = 1 and 7(11) = 0. It can be indicated as the regular one, 
since det(gik) = —L A = — 1 and can not vanish. Equation for /3 reduces to 



H 



(P'f [l+^ (A) 2 ] = (0') 2 +V(a) (0 2 )" + (A' 2 e-e) 2 + (A!^)' 



(73) 



When 77(A) — there is no real solutions of this equation, but such a possibility appears 
in the non-minimal case. We consider only one example of the exact regular models, it is 
characterized by 



= O , A 2 (u) = A 2 (0)e^ , A 3 (u) = A 3 (0)e~^ 
and is possible, when 77(A) < and 



0o |r/( A )| = 1 + ^(0)+ ^(0) 
The function (3(u) is arbitrary; we prefer to use the following periodic finite function: 

/3( M ) = ^/3 (max) (l-cos2A M ), 0(0) = 0, /3'(0) = 0. 



(74) 



(75) 



(76) 



The metric for this non-minimal model is regular: 

ds 2 =2dwdt> — jexp 2/3( max ) sin 2 Am (da; 2 ) 2 +exp — 2/3( max) sin 2 \u (dx 3 ) 2 } , (77) 
the potentials of the electromagnetic field and their derivatives are also periodic and regular. 

3.3.2 Non-minimal models with symmetric spacetime of the pp-wave type 

The key equation for the metric coefficients can be reduced in this case to the relation 



(A1 + A2) 1 + K77(A 



t/\2 



1+277( A) ) +2r/ (A) 00"+-l [(^e^)" + (j^e 



(78) 



Let us consider three exact solutions of this equation, which can be indicated as the non- 
minimal ones. 

(1) First exact solution 

Let us suppose that Ai = — A 2 = z/ 2 , and we again deal with the Petrov solution ( 164)) . This 
solution is admissible when 



|(0') 2 (l+277(A)J +2t/(a)00 // } cos 2 uu cosh 2 uu+ (A' 2 cosh vu) 2 + L4.3 cos 2 vu 



0. (79) 
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One of the solutions of this equation, which relates to the case V(A)=—\, is given by 



E 



4>(u)=(j)Q coshau , 0o= ± v^2 — , (80) 

a 

A2{u)=—^ — [lo smuu cos vu— vsinvucosuu] , (81) 

(ur— v l ) 

Az{u) = -- — rr [z/sinco'Msinhz/-u— a; cos cum cosh z/m] . (82) 

Here E , o~, u and v are arbitrary constant. When u 2 = u 2 , the expression for ^(m) should 
be replaced by .4.2 (w) = ^(2uu + sm2vu). Let us mention that in the paper [70], the case 
£ = 1/4 is also indicated as the special one (77(A) in our notations corresponds to £ in [70]). 

(2) Second exact solution 

In the case Ai = A2 = — /i 2 , an appropriate example of the model is characterized by 

0(10 = 0, L(u) = F(u) = G(u) = cosh/iw, L(0) = 1 , L'(0) = 0, (83) 

R u2u = R u3u = 7" = "Z 1 > ( 84 ) 

£■ 

^2(m)=-t-5 h~\ sin um cosh sinh /iw cos u«] , (85) 

£■ 

^4 3 (m)=— — — |/t sin uiu sinh /i-u— a; cos uu cosh , (86) 

{u +fl ) 

(f)(u) = 0o = const , 77( A ) < , kEI = 4/i 2 (\t] {k) \k4>1 - l) > . (87) 

Again, we deal with circularly polarized electromagnetic wave, whose amplitude E is con- 
nected with constant pseudoscalar field by a special manner (l8"Tj) . which is admissible at the 
non-minimal case only. 

(3) Third exact solution 

When Ai = u 2 , X 2 = — u 2 , there is a simple non- minimal modification of the solution (I6"8"|) . 
(|69|) with constant non-vanishing pseudoscalar field 

F(u) = cos uu , G(-u) = cosh /i-u , 0(w) = O , (88) 

" 2 - " 2 + 2(1 • (89 » 

with electromagnetic potentials, given by (170]) . ( ITT]) . 
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3.3.3 Special model 



When a pseudoscalar field is constant, i.e. <j){u) = 4>o, and electromagnetic field is absent, 
the equation for the gravitational field takes the form 



^+(/3') 2 cosh 2 2 7 + {if 



l+K?7(A)0o 



=0. 



(90) 



When 0o 



— 1/m7(a), the equations for gravity field are satisfied for arbitrary L(u), f3(u) 
and j(u), and we deal with a non-trivial special case. It corresponds to the case when the 
Ricci scalar disappears from the Lagrangian. 



3.3.4 Cheshire smile 

When the spacetime is the Minkowski one, and L = l,/3 = 0,7 = 0, the equation ( 1561) yields 

o= (<P'f +v(a } (<p 2 )" + {A? + (A) 2 , (9i) 

and this equation admits real solutions when the non-minimal parameter r/(A) is non- vanishing. 
The exotic situation when (pure) scalar field does not curve the spacetime, is called in [7T] 
as gravitational Cheshire smile. Here we deal with pseudoscalar field instead of pure scalar 
one; nevertheless, the physical sense of such a solution remains the same as in [71]. When 
the electromagnetic field is described by the circularly polarized wave with 

E E 
^4.2 (m) = — cos Xu , Az{u) = — sin Xu , (92) 
A A 

then the pseudoscalar field satisfies the equation 

(0f +?7(A)(0 2 ) // +£ 2 = O, (93) 
whose solution can be represented in quadratures as 

£2 



'C x <j> "( A > - 



1+27J (A) 

One of the explicit particular solutions to ( liMj) is 

1 E 

4> = 4>o cosh ku, 77 (A ) = --, <p G = y/2—. (95) 

Again we used the special value of the non-minimal parameter r]<x) = — | as in [70]. This 
solution is infinite at u — > oo. 

In table 2, we summarize and classify exact solutions obtained above for the model with 
pp-wave symmetry. 
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V(A) = 


»7(A) ^ 


4>(u)=4> 

00 7^0 


GW of Petrov's type (1M1) 


Regular periodic GEMW fTTl]) 
Regular aperiodic GEMW (J83]) 
GW with periodic EMW ([88]) 
Special solution (!90l) 


<f>'(u) ± 


Symmetric spacetime of GW type 
with axion field linear in time (1661) 


Axion- GEMW of Petrov's type flHD]) 
Axion-EMW Cheshire smile fl9T]) 


<f>{u) = 


Co-moving GEMW ([68]) 


Co-moving GEMW (J68]) 



Table 2. All the obtained solutions are classified with respect to the pseudoscalar field type: 
is a non- vanishing constant <fto; (ft depends on retarded time and (ft'(u) ^ 0; (ft vanishes (see 
horizontal lines); as well as, with respect to values of the guiding parameter: r/(A)=0 or ^ 0. 
We distinguish solutions of pure gravitational wave type (GW), co-moving gravito-electromagnetic 
waves (GEMW), as well as axion-GEMW and axion-EMW. The cases indicated as the solutions for 
symmetric spacetime relate to the model with covariantly constant Riemann tensor, the solutions 
of the Petrov type being the particular subcase of such spacetime with vanishing Ricci tensor. In 
parentheses we indicate the initial formula in the text related to the corresponding solution. The 
solutions with (ft(u) = are not new; they are introduced to complete the table. The solutions with 
r 7(A) = have, formally speaking, the same structure as the known minimal solutions for the true 
scalar field; we focus our attention to these solutions since they are also valid for the pseudoscalar 
field with arbitrary non- vanishing non-minimal coupling parameters q±, </3, Qi Q2, Q3, ^l, V2, 
7/3, thus belonging to the class of non-minimal solutions. 

4 Propagation of electromagnetic waves 
coupled to the background axion and 
gravitational pp-wave fields 

4.1 Evolutionary equations for the potentials 

Consider now test electromagnetic waves propagating in an arbitrary direction in the pure 
pp-wave background with 7(w) = 0. The test electromagnetic field satisfies the non-minimal 
evolutionary equations ( 120]) . f l2T]) with a vanishing current on the right-hand side, and the 
solutions, describing the axion field <ft(u) and pp-wave gravitational field, are assumed to 
be unperturbed. We deal below with pure gravitational waves; thus, the Ricci scalar and 
the Ricci tensor vanish (R = 0, R uu = 0), providing that the nonvanishing components of 
the Riemann tensor are of opposite signs (R 2 u2u = —R 3 U 3u)- Also, we require the potential 
four-vector to satisfy the Lorentz gauge condition 

21/ 

V k A k = d v A u +d u A v +g a0 d a Ap+— A v =0 . (96) 

Let us mention that the problem of propagation of electromagnetic waves coupled to axions 
in the geometric optics approximation is studied in the paper [72] on the basis of investigation 
of the Fresnel equation. Here we do not restrict ourselves by this approximation and consider 
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exact solutions of the equations of axion electrodynamics in analogy with solutions of the 
model of non-minimal optical activity discussed in [73] . 

Due to condition (|96l) the first equation from ( J201) (with i = u) accounting for (l2li) gives 
the following simple equation for the longitudinal potential A v = G v )Ai' 



VA, 



2L' 



f) A 



where T> is the differential operator: 



V = g mn d m d n = 2d u d v -L~ 2 (e^d 2 + e ^ d A . 



(97) 



(9£ 



Clearly, ( )9~7|) is an equation for A v only; it does not include other components of the potential 
four- vector. The exact solution of this equation takes the form (see [73] ) 

(99) 



A v = L-\u) B V {W), 
where B V (W) is an arbitrary function of the phase scalar W, given by 



Li r\jr>) 



+k v V+k2X 2 +k3X £ 



k\ J duL- 2 e- 2f} + k\ J duL- 2 e 2 ? 
with arbitrary constant values of Wo, k 2 , k 3 an d kv The wave four- vector Ki 

K . = ViW = _ 6 u_L ( g 22 k 2 + g 33 k 2^ + + ^ + ^ ? 



(100) 



(101) 



is the null four-vector, i.e. g im KiK m = and the longitudinal component of the electro- 
magnetic potential A v propagates with the speed of light in vacuum. A special solution 
B V (W) = has been motivated in |74j. This solution relates to the Landau gauge condition 
£L\Ai = 0, and we will use it below. Equations for the transversal components A 2 and A 3 
form a coupled system 



V + 2q 3 R 2 u2u d 2 -2(3'd v 



A 2 + e 



2fi 



<>'d v + 2Q^R 2 u2u d\ 



V + 2q 3 R\ 3u d 2 v + 2p'd v ] A 3 - [<j>'d v + 2Q 3 <f ) R 3 u3u d 2 
The replacement 

A % = ePB 2 , A 3 = e~?B 3 
simplifies these equations yielding 



^3 = 0, 

A 2 = 



V + 2q 3 R 2 u2u d 2 



V + 2q 3 R\ 3u d 2 



B 2 + 



B 3 



R + 2Q 3 <pR 2 u2u d 2 v 
i>'d v + 2Q 3( j ) R\ 3u d 2 v 



B 3 = 



Bo = 0. 



(102) 
(103) 

(104) 

(105) 
(106) 



Below we obtain and discuss exact solutions of these equations, suggesting that the functions 
<p(u), L(u) and /3(u) are presented by one of the exact solutions listed in the previous Section. 
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4.2 Birefringence induced by curvature in the case 
of vanishing pseudoscalar field ((f) = 0) 

In order to clarify a physical meaning of the parameter q 3 , let us consider the simplest model 
with = 0, i.e. the model without axion-photon interaction. As it was shown in [73], the 
solution of ([102]) and (TT031) at <p = are 



A v = 



A 2 = e^B (2) (W {2) ), 
1 



A 3 = e-% 3) (W {3) ), 



(107) 
(108) 



k2e- p B i2) (W i2) ) + k 3 e^B {3) (W (3) 
Here £>(2) and £>(3) are arbitrary functions of their arguments, W( 2 ) and W( 3 ), given by 

Wfa = W - q 3 k v £ duR 2 u2u (u) , W (3 ) =W- q 3 k v £ duR 3 u3u (u) , (109) 

with the phase W found in f llOOp . The numbers (2) and (3) indicate here the electromag- 
netic waves with polarization vectors directed along Ox 2 and Ox 3 axes, respectively. For 
estimations it is convenient to consider a particular case with harmonic functions 



B 



(2) 



B° {2) COS W(2) 



B 



(3) 



Bf 3) cos W {3) 



and parameter k v 



k,, 



W(EM) 

cy/2 



'1— COS I 



110) 



(111) 



explicitly expressed in terms of the electromagnetic wave frequency W(em) and angle 9 be- 
tween the directions of propagation of the gravitational wave (Ox 1 in our case) and of the 
electromagnetic one. 

Clearly, W(2) 7^ W(3), when the corresponding components of the curvature tensor do 
not coincide, R 2 u2u ^ R 3 u3u , and thus we deal with the birefringence effect induced by the 

(2) (3) 

curvature. The wave four-vectors Kr> = ViW (2) and K\ 0> = ViW {3) can be represented as 
follows: 



K 



(2) 



Ki - 5?q 3 k v R 



2 K 

u2u 5 ly 



(3) 



Ki - 6?q 3 k v R< 



u3u ' 



(112) 



where Ki is given by fllOip . The wave four- vectors k[ 2 ^ and Kf } are not longer null four- 
vectors, since 



-(3) 



g u K (2) K (2) = - 2 q 3 k 2 v R 2 u2u , fK^K^ = -2q 3 k 2 R 3 u3u . (113) 

The frequencies of the waves with polarization along Ox 2 and Ox 3 , respectively, (see, e.g., 
[75] for details), are 



CO 



(2) = 71 



^{g 22 k 2 2+ g 33 kt)-qa- 



u2u 



W( 3 ) 



v/2 



1 

2^2 



U3U 



(114) 
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These frequencies differ when R 2 u3u ^ R 3 u3u - The refraction indices n( 2 ) and 7i( 3 ) can be 
calculated as follows (see, e.g., [75] for details): 

\W(2)/ 



C 



2 




4) = 1 - — 9 ll Kr>Kr = 1 + 2g 3 — R\ 3u • (115) 

The difference An = n^ 2 )~ n {3) can t> e calculated for all model mentioned in the previous 
section, but we focus now on one example only. Let the test electromagnetic wave propagate 
in the direction — x 1 , and the solution for the gravitational field be of the Petrov type, i.e. 
R 2 U 2u = ~R 3 u3u = v<1 ■ Then one obtains the explicit expressions for the phase velocities: 

(2) (3) , 

^ = ^ = 1-^ 2 , ^ = ^ = 1W, k^^=, (116) 
c kc c kc v2 

and for the refraction indices: 

1 + q^v 2 1 — q^v 2 4g 3 z/ 2 

n (2) = \ 2 ' n (3) = T~i 2 ' An = i 2~T ' ( 117 ) 

1 — q 3 u 2 w 1 + q 3 u 2 1 — g 3 z/ 4 

This means, first, that the phase velocities differ from each other, one of them being less than 
the speed of light in vacuum c, while the second being bigger than c. Second, the non-minimal 
(curvature induced) difference between the refraction indices for the waves with orthogonal 
polarizations depends on the curvature v 2 and the non-minimal coupling parameter q 3 . Thus, 
in this illustration, we deal with birefringence in the absence of polarization rotation. 

The constraints on the non-minimal coupling parameter from solar system observations 
and from pulsar data analysis were studied in [75] (the parameter A in this paper relates 
to g 3 /4 in our case). From the timing observations of the binary pulsar PSR B1534+12 
signals, Prasanna and Mohanty gave a constraint on q 3 to be |g 3 | < 2.4 ■ 10 11 cm 2 . From 
more precise observations of the double pulsar system PSR J0737-3039 A/B, the Shapiro 
time delay (passing PSR J0737-3039 B of mass 1.25 Ms un and radius 10 km) is measured to 
be 6.21 ±0.33 /xs [77]. The theoretical estimate using general relativity is 6.153 ±0.26 /is. If 
we take an upper limit of 0.1 /is for the maximum contribution of the polarization-dependent 
correction, we have an upper bound |g 3 | < 2.4 • 10 10 cm 2 . 

4.3 Quasi-minimal model of optical activity (qs = 0, Q% = 0) 
induced by varying pseudoscalar field ((f) ^ (J)q) 

Now we consider the model in which the birefringence is absent, but the rotation of polar- 
ization is present. Let two non-minimal coefficients g 3 and Q 3 be vanishing; then the exact 
solutions of p>j). (I97jl . (TT02]) . (TT03|) are 

A 2 = e /3(M) | J B(^)cos^[0(M)-0(O)]±F(^)sin^[0( M )-0(O)]| , (118) 
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(119) 



(120) 



Here E'(W) and F(W) are arbitrary functions of one phase scalar W ( jlOOl) : thus, there is no 
birefringence. But now we deal with the effect of Faraday rotation, induced by the photon- 
axion interactions. This is a kind of optical activity induced by the axion field. The effect of 
the Faraday rotation is described by the rotation phase <j)(u)—<j)(0), so this effect disappears 
when the pseudoscalar field is constant. When the pseudoscalar field is a linear function 
of retarded time (see ( |66|) . (1671)). the Faraday rotation is uniform. Let us mention again 
that we indicate this model as quasi-minimal since we suppose that only two non-minimal 
parameter are vanishing, but other ones are non-vanishing. 

4.4 Non-minimal optical activity induced by photon-axion 

interaction in the case of a constant pseudoscalar field = (J)q 

When non-minimal coupling constants are non-vanishing, i.e. q 3 7^ and Q 3 7^ 0, we search 
for exact solutions to the system f l 1 5 f) . ( 11Q6|) as follows: 



Again W is given by ( llOOp . and four functions of the retarded time hi, h 2 , h 3 and hi satisfy 
the following system of linear ordinary first-order differential equations with a coefficient 
depending on retarded time u: 



Let the function <f>(u) be presented by a constant value <f>(u) = 0o- Introducing a new 
convenient variable r 



B 2 = hi(u)cosW + h 2 (u)smW, B 3 



h 3 (u) sin W + hi (u) cos W . 



(121) 




(122) 




(123) 



we reduce the system (11221) to the following system with constant coefficients: 



h i(r) = -q 3 h 2 + Q 3 4> h 



••1 • 



Q30o^i + <?3^4 , K(t) = -Q 3 4> h 2 - q 3 h 3 , 



(124) 
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where the dot denotes the derivative with respect to r. The characteristic equation for this 
system of equations with constant coefficients takes the form 

(A 2 + Q^ + g 3 2 ) 2 = 0, (125 ) 
and its solutions A = ±ifl are the double ones, where 

= \[ql + QWo ■ (126) 
Let us mention that (II 24ft in this case can be rewritten in terms of second-order equations: 

h a + Q 2 h a = 0, (127) 
where a = 1, 2, 3, 4. That is why the exact solution of the system (11241) 

hi{r) = fci(O) cos fir - — sinfir[g 3 h 2 (0) - Q 3 4> h 3 (0)] , 

h 2 (r) = h 2 (0) cos fir + - sin fir [q 3 hi(0) + Q 3 4>o h A (0)] , 
h(r) = h 3 (0) cos fir + ^sinfir[g 3 h 4 (0) - Q 3 (f> h^O)] , 

h 4 (r) = /i 4 (0)cosfir- ^sinfir[g 3 /> 3 (O) + Q 3 O h 2 (0)} , (128) 

contains the trigonometric functions cos fir and sin fir only Thus, the quantity fi plays 
a role of frequency of the curvature-induced rotation of the Faraday type. Let us stress 
that in this model the non-minimal coupling parameters Q 3 and g 3 play equivalent roles, 
predetermining the pure rotation of the polarization vector without damping. When the pp- 
wave spacetime is symmetric, i.e. the component R 2 u2u is constant, one obtains immediately 
from (I123p that r is linear in the retarded time, e.g. for solution ( 164")) one obtains r = v 2 k v u. 

4.4.1 Constraints on the non-minimal coupling parameters from CMB polar- 
ization observations 

The parameter Q 3 enters the key formulas in the product Q 3 0o and thus relates to the non- 
minimal axion contribution to the effect of polarization rotation. Let us assume that g 3 = 
and the following simple initial conditions take place 

MO) = /is(0) = /i4(0) = 0, h 2 (0)=B ->• B 2 = B smW, B 3 = , (129) 

i.e. initially the electromagnetic wave was polarized along the Ox 2 axis only. Then we obtain 
from ( 1128]) 

A 2 = B e p cos (Q 3 <Pot) sin W , A 3 = -B Q e~ p sin (Q 3 <P t) cos W . (130) 
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Clearly, when Q 3 (f)o 7^ 0, the orthogonal component of the electromagnetic potential A 3 
appears; thus, the polarization rotation can exist when the parameter Q 3 (po is non- vanishing. 
The polarization ellipticity can be described in this case by the formula 



(131) 



cos 2 (Qs^ot) ' sin 2 (Q 3 (j)oT) 
Thus, the polarization rotation rate is predetermined by the function Q 3 (fio T , an d the quantity 
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can be considered as an equivalent of the gradient of the pseudoscalar field for the estimates 
of the effect. The frequency dependence is different from the Faraday rotation in an ionized 
medium. From f II 3 2 [) . the order of magnitude of [V0] e ff e ctive from a microwave of frequency 
100 GHz is 

[V0] cffcct i vc ~ 4.7 Q 3 (p 



'(GW) 
X 2 

A (GW) 



cm 



hence, 



Let H, 



(GW) 



|A0 effec tive| - 4.7 cm 1 Q 3 4>o 
10~ 15 £, where £ < 1; then for 
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(136) 



Constraints on cosmic polarization rotation angle (p come from three kinds of observations: 
(i) polarization observations of radio galaxies; (ii) optical/UV polarization observations of 
radio galaxies; (iii) CMB polarization observations. Radio observations put a limit of Ap < 
0.17 — 1.0 rad over cosmological distance from various types of analysis (see, e.g., [21] for 
a review). Optical/UV observations put a limit of Aip < 0.17 rad [7H1 [79]. The QUaD 
observations [SU] on CMB polarization give the most stringent constraint on A(p 



Aip = 0.0096 ± 0.0143 ± 0.008 rad . 



(137) 



For gravitational waves with a wavelength of the order of the Hubble distance, this is the 
constraint, A<£>hd; for shorter wavelength, A, we have to use the estimate 



Apx ~ Ap 



A 



HD 



:i38) 



Hubble Distance 

Constraints on polarization rotation 
data are under study. 
2tp, we have the con- 



In the general case, Ap depends on directions (22 

could also come from solar system observations: the analysis of these data are under study. 
Combining the estimate (I136p with f ll38p and using the relation 
straint Q 3 (f) < 2 • lO^^cm 2 . 
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5 Conclusions 



1. We formulated a new non-minimal ten-parameter Einstein-Maxwell-axion model, i.e. on 
the basis of the Lagrangian approach, we derived a non-minimally extended self-consistent 
(coupled) system of equations for electromagnetic (see equations (l20l) . fr2T|) . (!22|) ). pseudoscalar 
(see equation (1231) ) and gravitational (see (I24l) - (l3~71) ) fields. The Lagrangian is linear in the 
curvature tensor and its contractions, and is quadratic in the Maxwell tensor, i.e. we deal 
with one of the versions of non-minimal linear axion electrodynamics. 

2. The constitutive tensor (see (HUl) with (fl4l) - (Tl8l) ). associated with this model, is manifestly 
symmetric with respect to transposition of pairs of indices, i.e. this model excludes skewons 
[27] and describes non-minimal interaction of photons and axions only. An extension of the 
model accounting for skewon-type interactions will be considered separately. 

3. Decomposition of the constitutive tensor demonstrates explicitly that the dielectric 
permittivity, magnetic impermeability tensors and the tensor of cross-effects (or tensor of 
magneto-electric coefficients, in other words) acquire non-minimal contributions, including 
terms proportional to the pseudoscalar (axion) field. Since the tensor of non-minimal cross- 
effects is non- vanishing, one may conclude that, generally, the curvature-induced effect of 
optical activity is expected, even if the pseudoscalar (axion) field is constant. 

4. We applied the non-minimal Einstein-Maxwell-axion model to the case when the space- 
time has the pp-wave symmetry and all the physical fields inherit this symmetry. The 
reduced system of equations contains in this case only one key-equation (1551) . We presented 
nine examples of exact solutions for the non-minimally coupled electromagnetic, axion and 
gravitational fields with pp-wave symmetry. One can emphasize three new solutions among 
them. First, there is a new regular solution (see section 3.3.1.), which can appear in the 
non-minimal model only; the term regularity means that gravitational, electromagnetic and 
pseudoscalar fields have no singularities and are described by the functions finite everywhere. 
Second, there exist exact solutions describing circularly polarized electromagnetic waves cou- 
pled non-minimally to the axion field and plane gravitational waves (see (I80l) - (l82l) ). Third, 
the solution which is known as Cheshire smile is also admissible in this model (see section 
3.3.4). 

5. Discussing the optical activity induced by photon-axion interaction, we solved exactly the 
problem of propagation of test electromagnetic waves coupled non-minimally to the pseu- 
doscalar and gravitational fields in the pp-wave background. According to the classification 
of the models of photon-axion interaction given in the paper [72] , we deal here with the case 
when the four-gradient of the pseudoscalar field is a null four-vector, i.e. Vfc0V fe = 0. We 
show explicitly that the non-minimal coupling of the photon and axions with gravitational 
field generally leads to the birefringence effect. We discussed explicit exact solutions which 
describe the effect of optical activity. 

6. Non-minimal Einstein-Maxwell-axion model can give the cosmological solutions of the 
FLRW and Bianchi-I type, as well as static spherically symmetric solutions of the Reissner- 
Nordstrom type. We intend to consider the corresponding exact solutions in the near future. 



24 



Acknowledgments 

We thank the National Natural Science Foundation of China (Grant No. 10875171) and 
Russian Foundation for Basic Research (Grants No. 08-02-00325-a and 09-05-99015) for 
support. A.B. is grateful to colleagues from the Center for Gravitation and Cosmology of 
Purple Mountain Observatory of Chinese Academy of Science for hospitality. 

References 

Kallosh R, Linde A, Linde D and Susskind L 1995 Gravity and global symmetries Phys. 



Rev. D52 912 (arXiv:hep-th/9502069) 



Peccei R D 2008 The Strong CP problem and axions Lect. Notes Phys. 741 3 
( |arXiv:hep-ph/0607268[ ) 

Raffelt G G 1990 Astrophysical methods to constrain axions and other novel particle 
phenomena Phys. Rep. 198 1-113 

Turner M S 1990 Windows on the axion Phys. Rep. 197 67-97 

Battesti R et al 2008 Axion searches in the past, at present, and in the near future Lect. 
Notes Phys. 741 199 flarXiv:0705.0615D 

Ni W-T 1973 Nonmetric theory of gravity preprint (Montana, USA: Montana State 



University, Bozeman) (http://cgc.pmo.ac.cn) 



Ni W-T 1974 Weak equivalence principles and gravitational coupling Bull. Am. Phys. 
Soc. 19 655 

Ni W-T 1977 Equivalence principles and electromagnetism Phys. Rev. Lett. 38 301 

Peccei R D and Quinn H R 1977 CP conservation in the presence of pseudoparticles 
Phys. Rev. Lett. 38 1440 

Weinberg S 1978 A new light boson? Phys. Rev. Lett. 40 223 

Wilczek F 1978 Problem of strong P and T invariance in the presence of instantons 
Phys. Rev. Lett. 40 279 

Wilczek F 1987 Two applications of axion electrodynamics Phys. Rev. Lett. 58 1799 

Breton N 2000 Exact solutions in Einstein-Maxwell-dilaton-axion theory Recent devel- 
opments in gravitation ed J Ibanez (Espana: Univ. Pais Vasco Proceedings of ERE 99) 
p 179 

[14] Sur S, Das S and SenGupta S 2005 Charged black holes in generalized dilaton-axion 
gravity JEEP 0510 064 (arXiv: hep-th/050815) 



25 



Bassett B A and Kunz M 2004 Cosmic acceleration vs axion-photon mixing Astrophys. 



J. 607 661 (arXiv: astro-ph/0311495) 



Halilsoy M and Sakalli I 2003 Collision of electromagnetic shock waves coupled with 



axion waves: An Example Class. Quantum Grav. 20 1417 (arXiv: gr-qc/0303102) 



Nojiri S, Odintsov S D, Ogushi S, Sugamoto A and Yamamoto M 1999 Axion- 
dilatonic conformal anomaly from Ads/CFT correspondence Phys.Lett. B 465 128 
( |arXiv:hep-th/9908066D 



Lopez L A and Breton N 2007 Asymptotic singular behaviour of inhomogeneous cos- 
mologies in Einstein-Maxwell-dilaton-axion theory Gen. Relat. Grav. 39 153 (arXiv: 
|gr-qc/0608125l ) 

Carroll S M, Field G B and Jackiw R 1990 Limits on a Lorentz- and parity-violating 
modification of electrodynamics Phys. Rev. D 41 1231 

Carroll S M and Field G B 1991 Einstein equivalence principle and the polarization of 
radio galaxies Phys. Rev. D 43 3789 

Ni W-T 2008 From equivalence principles to cosmology: cosmic polarization rotation, 
CMB observation, neutrino number asymmetry, Lorentz invariance and CPT Prog. 
Theor. Phys. Suppl. 172 49 flarXiv:0712.4082p 

Ni W-T 2009 Cosmic polarization rotation, cosmological models, and the detectability 
of primordial gravitational waves Int. J. Mod. Phys. A 18 & 19 3493 (jar Xiv :0903. 0756) 

Zavattini E et al (PVLAS Collaboration) 2006 Experimental observation of optical 
rotation generated in vacuum by a magnetic field Phys. Rev. Lett. 96 110406 

Chen S-J, Mei H-H and Ni W-T 2007 Q & A experiment to search for vacuum dichroism, 
pseudoscalar-photon interaction and millicharged fermions Mod. Phys. Lett. A 22 2815 
( |arXiv:hep-ex/0611050p 

Zavattini E et al 2008 New PVLAS results and limits on magnetically induced optical 
rotation and ellipticity in vacuum Phys. Rev. D 77 032006 (arXi v:0706. 3419) 

Battesti R et al 2008 The BMV experiment: a novel apparatus to study the propagation 



of light in a transverse magnetic field Eur. Phys. J. D 46 323 ( arXiv:07 10.1703) 



Hehl F W and Obukhov Yu N 2003 Foundations of classical electrodynamics: Charge, 
flux, and metric (Boston: Birkhauser) 

Hehl F W and Obukhov Yu N 2005 Linear media in classical electrodynamics and the 



Post constraint Phys. Lett. A 334 249 (arXiv:physics/0411038) 



Hehl F W and Obukhov Yu N 2005 Measuring a piecewise constant axion field in 
classical electrodynamics Phys. Lett. A 341 357 ( jarXiv: physics / 0504 1 72 ) 



26 



Faraoni V, Gunzig E and Nardone P 1999 Conformal transformations in classical gravi- 



tational theories and in cosmology Fundam. Cosm. Phys. 20 121 (arXiv:gr-qc/9811047) 



Amendola L 1993 Cosmology with non-minimal derivative couplings Phys. Lett. B 301 
175 

Capozziello S and Lambiase G 1999 Non-minimal derivative coupling and the recovering 
of cosmological constant Gen. Rel. Grav. 31 1005 ( |arXiv:gr-qc/9 901051) 

Capozziello S, Lambiase G and Schmidt H-J 2000 Non-minimal derivative couplings and 
inflation in generalized theories of gravity Annalen Phys. 9 39 ( |arXiv:gr-qc/ 9906051) 

Prasanna A R 1971 A new invariant for electromagnetic fields in curved spacetime Phys. 
Lett. A 37 331 

Hehl F W and Obukhov Yu N 2001 How does the electromagnetic field couple to gravity, 
in particular to metric, nonmetricity, torsion, and curvature? Lect. Notes Phys. 562 
479 ( arXiv:gr-qc/0001010[ ) 

Novello M and Salim J M 1979 Nonlinear photons in the universe Phys. Rev. D 20 377 

Drummond I T and Hathrell S J 1980 QED vacuum polarization in a background 
gravitational field and its effect on the velocity of photons Phys. Rev. D 22 343 

Accioly A J, Vaidya A N and Som M M 1983 Non-minimal coupling and Bianchi type-I 
cosmologies Phys. Rev. D 28 1853 

Goenner H F M 1984 Theories of gravitation with non-minimal coupling of matter and 
the gravitational field Found. Phys. 14 865 

Turner M S and Widrow L M 1988 Inflation-produced, large-scale magnetic fields Phys. 
Rev. D 37 2743 

Mohanty S and Prasanna A R 1998 Photon propagation in Einstein and higher derivative 
gravity Nucl. Phys. B 526 501 QarXiv:gr-qc/9804017[ ) 

Teyssandier P 2004 Variation of the speed of light due to non-minimal coupling between 
electromagnetism and gravity Annales Fond. Broglie 29 173 (arXiv:gr-qc/0303081) 



Lafrance R and Myers R C 1995 Gravity's rainbow: Limits for the applicability of the 
equivalence principle Phys. Rev. D 51 2584 (jar Xiv:hep-th/ 941 1018) 

Solanki S K et al 2004 Solar constraints on new couplings between electromagnetism 
and gravity Phys. Rev. D 69 062001 ( arXiv:gr-qc/0402055[ ) 

Kostelecky A and Mewes M 2002 Signals for Lorentz violation in electrodynamics Phys. 
Rev. D 66 056005 flarXiv:hep-ph/0205211[ ) 



27 



Bamba K and Odintsov S D 2008 Inflation and late-time cosmic acceleration in non- 
minimal Maxwell- .F(.R) gravity and the generation of large-scale magnetic fields JCAP 
0804 024 (larXiv:0801 .095411 

Balakin A B and Lemos J P S 2005 Non-minimal coupling for the gravitational and 
electromagnetic fields: A general system of equations Class. Quant. Grav. 22 1867 
flarXiv:gr-qc/0503076D 

Balakin A B, Bochkarev V V and Lemos J P S 2008 Non-minimal coupling for the 
gravitational and electromagnetic fields: black hole solutions and solitons Phys. Rev. D 
77 084013 ( arXiv:0712.4"066l 

Rubilar G F, Obukhov Yu N and Hehl F W 2003 Torsion non- minimally cou- 
pled to the electromagnetic field and birefringence Class. Quantum Grav. 20 L185 
( |arXiv:gr-qc/0305049p 

Itin Ya and Hehl F W 2003 Maxwell's field coupled non-minimally to quadratic tor- 
sion: Induced axion field and birefringence of the vacuum Phys. Rev. D 68 127701 
darXiv:gr-qc/0307063[ ) 

Balakin A B and Zayats A E 2007 Non-minimal Wu-Yang monopole Phys. Lett. B 644 
294 ( |arXiv:gr-qc/0612019[ ) 

Balakin A B, Sushkov S V and Zayats A E 2007 Non-minimal Wu-Yang wormhole Phys. 
Rev. D 75 084042 flarXiv:0704. 12241) 

Balakin A B, Dehnen H and Zayats A E 2007 Non-minimal Einstein- Yang-Mills-Higgs 
theory: Associated, color and color-acoustic metrics for the Wu-Yang monopole model 
Phys. Rev. D 76 124011 f!arXiv:0710.5070ft 

Balakin A B, Dehnen H and Zayats A E 2008 Non-minimal pp-wave Einstein- Yang- 
Mills-Higgs model: color cross-effects induced by curvature Gen. Relat. Grav. 40 2493 
flarXiv:0803.3442p 

Bamba K, Nojiri S and Odintsov S D 2008 Inflationary cosmology and the late-time 
accelerated expansion of the universe in non-minimal Yang-Mills- _F(i?) gravity and non- 
minimal vector-F(i?) gravity Phys. Rev. D77 123532 flarXiv: 0803 .3384ft 

Balakin A B, Dehnen H and Zayats A E 2008 Effective metrics in the non-minimal 
Einstein- Yang-Mills-Higgs theory Annals of Physics 323 2183 flarXiv:0804.2196D 

Eringen A C and Maugin G A 1989 Electrodynamics of continua (New York: Springer- 
Verlag) 

Landau L D, Lifshitz E M and Pitaevskii L P 1996 Electrodynamics of continuous media 
(Oxford: Butterworth Heinemann) 

O'Dell T H 1970 The electrodynamics of magneto- electric media (Amsterdam: North- 
Holland) 



28 



Ni W-T 1983 Equivalence principles, their empirical foundations, and the role of preci- 
sion experiments to test them Precision measurement and gravity experiment ed W-T 
Ni (Taipei, Republic of China: National Tsing Hua Univ.) p 491 

Ni W-T 1984 Equivalence principles and precision experiments Precision measurement 
and fundamental constants II eds B N Taylor and W D Phillips (Washington D C, 
USA: National Bureau of Standarts special publication No 617 US Government printing 
Office) p 647 

Ni W-T 1984 Timing observations of the pulsar propagations in the galactic gravita- 
tional field as precision tests of the Einstein equivalence principle Proceedings of the 
Second Asian-Pacific Regional Meeting of the International astronomical Union eds B 
Hidayat and M W Feast (Jakarta, Indonesia: Tira Pustaka) p 441 

Lammerzahl C and Hehl F W 2004 Riemannian light cone from vanishing birefringence 
in premetric vacuum electrodynamics Phys. Rev. D 70 105022 ( |arXiv:gr-qc/ 0409072) 

Misner C W, Thorne K S and Wheeler J A 1973 Gravitation (San Francisco: Freeman) 

Balakin A B and Lemos J P S 2001 Singular behaviour of electric and magnetic fields in 
dielectric media in a non-linear gravitational wave background Class. Quantum Grav. 
18 941 ( arXiv:gr-qc/0101123D 

Bondi H, Pirani F A E and Robinson I 1959 Gravitational waves in general relativity, 
III, Exact plane waves Proc. Roy Soc. A 251 519 

Petrov A Z 1969 Einstein spaces (Oxford UK: Pergamon Press) 

Synge J L 1971 Relativity: the general theory (Amsterdam: North-Holland) 

Stephani H, Kramer D, MacCallum M, Hoenselaers C and Herlt E 2003 Exact solutions 
of Einstein's field equations (Cambridge, UK: University Press) 

Ayon-Beato E and Hossaine M 2005 Scalar fields non-minimally coupled to pp waves 
Phys. Rev. D 71 084004 flarXiv:hep-th/0501040D 

Ayon-Beato E, Martinez C, Troncoso R and Zanelli J 2005 Gravitational Cheshire effect: 
Non-minimally coupled scalar fields may not curve spacetime Phys. Rev. D 71, 104037 



(arXiv:hep-th/0505086) 



Itin Y 2008 Wave propagation in axion electrodynamics Gen. Relat. Grav. 40 1219 
f!arXiv:0706.2991h 

Balakin A B and Lemos J P S 2002 Optical activity induced by curvature in a gravita- 
tional pp-wave background Class. Quantum Grav. 19 4897 (arXiv:gr-qc/0209062) 

Balakin A B 1997 Gravitational radiation and birefringence induced by curvature Class. 
Quantum Grav. 14 2881 



29 



Balakin A B, Kerner R and Lemos J P S 2001 Cherenkov radiation in a gravitational 



wave background Class. Quantum Grav. 18 2217 (arXiv:gr-qc/0012 101[ ) 



Prasanna A R and Mohanty S 2003 Constraints on non-minimally coupled curved 
space electrodynamics from astrophysical observations Class. Quant. Grav. 20 3023 
( |arXiv:gr-qc/0306021[ ) 



Kramer M et al 2006 Tests of general relativity from timing the double pulsar Science 
314 97 (arXiv:astro-ph/060941) 

Cimatti A, di Serego Alighieri S, Field G B and Fosbury R A E 1994 Stellar and 
scattered light in a radio galaxy at z = 2.63 Astrophys. J. 422 562 

di Serego Alighieri S, Field G B and Cimatti A 1995 The physics of the interstellar 
medium and intergalactic medium ASP Conference Series 80 276 

Pryke C et al 2009 Second and third season QUaD CMB temperature and polarization 
power spectra Astrophys. J 692 1247 ( 1arXiv:0805. 1944ft 

Wu E-Y-S et al 2009 Parity violation constraints using cosmic microwave background 
polarization spectra from 2006 and 2007 observations by the QUaD polarimeter Phys. 
Rev. Lett. 102 161302 flarXiv:0811.0618p 



30 



